Motion-generated quantum gates and entanglement resonance 
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The motion of a cavity in spacetime generates quantum gates via entanglement between cavity 
modes. While entanglement produced by uniform acceleration is small, we show analytically that 
two-mode quantum gates with an observable amount of entanglement can be generated by periodic 
repetition of any non-uniform trajectory segment. A resonance occurs when the frequency associated 
with the duration of the repeated segment equals half the frequency sum of any two oddly separated 
cavity modes, leading to an entanglement accumulation that grows linearly in the number of segment 
repetitions. 



Introduction. — Understanding how motion and grav- 
ity affect entanglement is a key feature in the implemen- 
tation of new quantum information technologies, includ- 
ing quantum cryptography and teleportation, in space- 
based scenarios that are currently under investigation [1] . 
Quantifying entanglement in situations where motion 
or gravitation have a significant role can also provide 
guidance for theories about the microscopic structure of 
spacetime, via the Hawking-Unruh effect and its connec- 
tions to thermodynamics and statistical mechanics [3] . 
Recent work in relativistic quantum information has 
shown that non-uniform motion of a cavity creates en- 
tanglement between the cavity field modes [U [5] . In this 
Letter we show that by repeating any trajectory segment 
periodically it is possible to generate a resonance cre- 
ating higher amounts of entanglement for small cavities 
and large accelerations or large cavities and small acceler- 
ations. Non-uniform motion induces transformations on 
the states which are equivalent to two-mode squeezing 
operations upto local unitaries. These transformations 
play the role of two-mode quantum gates in continuous 
variables systems 0. 

Finding suitable ways to store and process information 
in a quantum and relativistic setting is a main goal in the 
field of relativistic quantum information. Moving cavities 
are good candidates to store information [7 - 9 since con- 
fined fields can be realized experimentally and observers 
can directly access their states by means of local opera- 
tions. When a cavity is accelerated for a finite time the 
cavity modes are affected by the motion. A mismatch 
between the vacua at different times gives rise to the cre- 
ation of particles which populate the modes [TO]. The 
initial and final modes are related through Bogoliubov 
transformations which mix all frequency modes |3J [7J 
generating entanglement [H [5]. Quantum correlations 
are created between all modes; however, higher degrees 
of entanglement are produced between oddly separated 
modes. As a consequence, maximally entangled states of 
modes confined in two different cavities degraded when 
the cavities undergo non- uniform motion [7J. 

We show that two-mode squeezed states can be pro- 



duced in a single cavity by the periodic repetition of 
any trajectory segment. Entanglement resonances occur 
when the frequency associated to the segment travel time 
is equal to the sum of two oddly separated cavity mode 
frequencies. We show analytically that for any couple 
of oddly separated modes it is possible to find a seg- 
ment travel time where the entanglement between such 
modes increases linearly with the number of repetitions. 
These resonances appear independently of the details of 
the trajectory though the amount of entanglement gen- 
erated does depend on the trajectory itself. 

As a concrete example we present a travel scenario 
which allows for simple analytical expressions. The travel 
segment which we repeat starts with a period of accel- 
eration followed by inertial coasting. The cavity then 
accelerates in the same or opposite direction for a sec- 
ond period of time and finally undergoes a second period 
of inertial coasting. We find the conditions for the reso- 
nances to occur and investigate how the magnitude and 
the direction of the accelerations affect the final entan- 
glement generated. In the special case where there is 
no inertial coasting and the accelerations alternate in di- 
rection, our sample trajectory reduces to the oscillatory 
motion that is often considered in the dynamical Casimir 
effect literature 

Setup. — We consider a real scalar field </> of mass m 
contained within a cavity in (1 + l)-dimensional space- 
time, with Dirichlet boundary conditions at the walls. 
The massless field can be treated as a special case of our 
study and the effect of additional transverse dimensions 
can be included as a positive contribution to m. The cav- 
ity follows a worldtube which is composed of periods of 
inertial and uniformly accelerated motion. We will start 
by describing the field within the cavity during such pe- 
riods as seen by a co-moving observer. 

During periods of inertial motion, we use Minkwoski 
coordinates (t, x) to describe our system. The cav- 
ity walls are placed at x = xa and x = xb where 
< xa < xb, so that L = xb — xa is the length of 
the cavity. The positive frequency mode functions with 
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respect to the time translation Killing vector <9 t are: 



<P k (t,x) = 



1 



sin 



kn 



(x - x A ) 



—iu> k t 



(1) 



where k £ N (we set c = h = 1). 

We employ Rindlcr coordinates (r, x) to describe the 
field during periods of uniformly accelerated motion. 
Rindlcr coordinates are related to Minkowski coordi- 
nates through the transformations x — six 2, — $ and 
t = (1/a) atanh(i/x), where a is a positive constant. An 
uniformly accelerated observer comoving with the cavity 
follows the boost Killing vector field d v = xdt+td x which 
in Rindler coordinates takes the form d r] — (l/a)d T . The 
cavity walls are placed at \ = xa and x = x b and the 
proper time and acceleration at the centre of the cavity 
are given by r and a — 2/ (xb+xa), respectively. The so- 
lutions to the Klein-Gordon equation that are of positive 
frequency with respect to d n are 

<l>k ( T , X) =N k [l l n k /a{mx)I^ l n k /a{ , mxA) 

-I-in k /a(mx)I t n k /a(mx A )]e- mkT , (2) 

where I v are modified Bessel functions of the first 
kind [12] and N k are normalisation constants. In the 
massless case the mode functions and the frequencies re- 
duce to simple expressions (see [7]). 

The quantized field operators are given by <fi M = 
Efcfak'a* + h.c), and ^ = EM A k + h.c.) dur- 
ing periods of inertial and accelerated motion, respec- 
tively. The Minkowski and Rindler annihilation and 
creation operators obey the standard commutation rela- 
tions [afc,a|] = 5^ and = 5m- The Minkowski 
frequencies, io k > 0, are determined by the Dirichlet 
boundary condition <p^(t, xb) — and given by ([lj, 
while Rindler frequencies, £lk > 0, are determined by 
cf>%(T,x B ) = 0. 

We will work in the covariant matrix formalism which 
is applicable to systems consisting of a discrete number 
of bosonic modes as long as the analysis is restricted to 
gaussian states. In this framework the state of the sys- 
tem is entirely described by its first and second moments 
[T31 ITi] . The evolution of the state is given by a simi- 
larity transformation S Tp AS where S is the symplectic 
representation of the evolution and A is the covariance 
matrix encoding all information pertinent to the state. 
Tp denotes matrix transposition. 

Travel scenario techniques. — Changes from inertial to 
accelerated motion and vice versa are implemented by 
the action of Bogolubov transformations. Consider that 
at t = a cavity initially at rest begins to accelerate. The 
inertial and accelerated cavity modes are related by the 
Bogolubov transformation </>f = J2 n { a kn(t>n + Pkn4>n*) > 
where the star denotes complex conjugation and the Bo- 
golubov coefficients a k i and can be evaluated by tak- 
ing Klein-Gordon inner products of the two sets of modes 



at t = [Till IHJ [IB]. In the low acceleration regime 
we can expand the coefficients in a Maclaurin series as 
a = l+ha^+0(h 2 ) and f3 = h(i^+0(h 2 ) where h = aL 
is a small dimensionless expansion parameter [7, 8j. The 
coefficients can be computed using uniform asymptotic 
expansions of Bessel functions [T^l [T7] . 

During the periods of inertial or accelerated motion 
the modes undergo free evolution which induces phase 
rotations on the state of the form U(t) — © i=1 R(0i) 
where N is the number of bosonic modes considered and 
R is the rotation matrix. 



R{6 k ) := 



cos u k sm t>k 
— sin 8k cos 9 k 



(3) 



The angles are given by 9 k = uikt during coasting peri- 
ods and 9 k = SlfcT during acceleration. The Rindler and 
Minkowski frequencies coincide to first order (flk(h) 
Lu k & as h — > 0). 

By composing these basic transformations we con- 
struct our trajectores. The field modes of a cavity ini- 
tially at rest and the modes after any travel scenario are 
related through general Bogolubov transformations. In 
the covariant matrix formalism the transformations are 
represented the symplectic matrix 
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where the 2x2 matrices Skw take the form 



Skk' 



Bkk') 3(^4fcfc' + Bkk' 



kk' ^kk' ) ^y J1 -kk' t J-'kk' ) \ /r\ 

%(A kk , - B kk .) $t{A kk , + B kk .) [0) 



where A kk > and B kk > are Bogolubov coefficients associ- 
ated with the whole trajectory [5]. In the case of a cav- 
ity initially at rest that begins to uniformly accelerate 
at t = the symplectic matrix, which we denote by V, 
corresponds to A kk > = a kk > and B kk > = hk- ■ When 
the Bogolubov coefficients f3 kk i are non- vanishing there is 
particle creation and the cavity modes become entangled 
according to a co-moving observer. The basic building 
block trajectory which corresponds to inertial-uniformly 
accelerated-inertial motion, is implemented by the action 
of the symplectic matrix = V _1 (/i)C/(r)V(/i). The en- 
tanglement generated between the cavity modes after a 
single basic building block trajectory has been analyzed 
in [5] and found to be very small. 

More general travel scenarios are constructed by com- 
posing basic building block trajectories. Here we are in- 
terested in constructing trajectory segments which will 
then be repeated to generate resonances. A segment con- 
tains any number of different basic building blocks and 
could be for example, a return voyage to Alpha-Ccntauri 
[7]. We consider the cavity to be initially in the vacuum 
state which satisfies a k \0) = Vk and is represented by 
the identity matrix A = I. This choice is convenient since 
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the state is gaussian and separable (not entangled) in the 
mode number degree of freedom k. In this case, it is pos- 
sible to calculate the reduced state of two field modes k 
and k' after the cavity has followed a general trajectory 
segment as long as the product of the acceleration and 
the cavity length remains small. We assume the different 
accelerations involved in the segment to be of the form 
di = Sja where a is the largest acceleration and < 1 
a real number. The Bogolubov coefficients can be ex- 
panded to first order in h = aL as A^k 1 — GkSkk' + A k x k \ 



and B 



kk' 



^kk' wnere the superscript (1) denotes the 



first order correction and Gi 



are the phases ac- 
quired by the state during periods of free evolution. T 
denotes the total proper time of the segment. When the 
paramater h is small and the initial state is the vacuum, 
the reduced state of two modes after a segment trajectory 
is given by <j\ = S kk ,S k k' and by a N = (S kk ,) Tp S kk , af- 
ter the segment has been repeated N times. The matrix 
Skk' takes the form 



Skk' — 



Skk Skk' 
Sk'k Sk'k' 



(6) 



Entanglement resonance and two-mode quantum 
gates — Our aim is to quantify the entanglement of the 
reduced state crjv- We find that, to first order in h, the 
final reduced state of two modes is separable for evenly 
separated modes but entangled and pure for oddly sep- 
arated modes since Tr(cri) = 1 + (D(h 2 ) [5]. The mixed- 
ness between k, k! and all other modes contributes only 
at second order 0(h 2 ) [5]. The entanglement of a bi- 
partite system in a pure state is quantified by the von 
Neumann entropy; however, this measure is not suitable 
in our perturbative regime because the first order con- 
tribution in the expansion of the entropy is of the form 
Svn ~ h In h which cannot appear as a term in a power 
series expansion. Fortunately, lower bounds to the en- 
tanglement generated in the system can be found by cal- 
culating the smallest positive symplectic eigenvalue vn 
of the partial transposed state &n = Po'nP Q3] where 
P = diag(l, 1, 1, —1). This symplectic eigenvalue is found 
by diagonalizing the matrix M — iQa^ where f2 is the 
antisymmetric symplectic form which, written in terms 
of Pauli matrices a y , takes the form Q = —Oy®o y . One 
can show that Vjq has a well defined power series expan- 
sion in h. There is a class of entanglement measures, 
including the logarithmic negativity, based on the posi- 
tive partial transpose criterion [18j which are monotonic 
functions of vm and therefore, also have a suitable power 
series expansion. 

The perturbative expansion of the symplectic eigen- 
value is given by Vjsr = 1 — ^jy where < vjy oc h is 
the first order correction. We find that when the com- 
mutator [S kk ,, Skk'] vanishes the partial transposed state 
after N segment repetitions is equal to the N-th power 
of the partial transposed state after a single repetition of 

~N 



is possible to show that to first order, when the commu- 
tator vanishes, the partial transpose state grows linearly 
with the number of repetitions 

After diagonalizing d\ and ctjv we find that the correc- 
tions to the eigenvalues satisfy the relation i>jy = Nv^. 
Any entanglement measure E(un) that is a function of 
the symplectic eigenvalue Djf will have a power expan- 
sion with pjy as a leading term. In particular, when 
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entanglement grows linearly 
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and therefore [S kk , , Skk'] = is a resonance condi- 
tion. In particular, the logarithmic negativity at reso- 
nance is given by Ejs/ = Ni>i. Note that the condition 
Ni>P <C 1 must be satisfied to remain within the pertur- 
bative regime [5]. The state we consider is pure, bipartite 
and symmetric. It is well known in the covariant matrix 
formalism that in this case the state is equivalent to a two 
mode squeezed state upto local operations which do not 
affect the entanglement 14J. Operations that produce 
two-mode squeezed states are used as two-mode quan- 
tum gates in continuos variable systems. Using this fact 



one can show that v^p ~ Bw where B kk , is the first 
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order correction to the beta coefficients B k k' of the seg- 
ment [S]. Therefore, a lower bound to the entanglement 
generated at resonance after N segment repetitions will 
be given by the logarithmic negativity Ejv = NB kk \ . We 
find that the commutator for an arbitrary segment to 
first order in h, is 



qTp o 
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where 



C n = 
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(cos(w fe T) - coB(wfc/T))K(B^) 
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C22 



(8) 



C12 



+ (sin(w fc T) +sm{u k ,T))%{B ( k l h 



(- cos(w fc T) + cos(w fc /T))9(B$) 



+ (sinKT) + sra(wvT))K(B$) 



(9) 



Since the Minkowski and Rindler frequencies coincide 
to order ft, 1 , the coefficients CV, depend only on the 
Minkowski frequencies cj/.. We find that the commutator 
vanishes when the following equation is satisfied 



B$[G£-Gv]=0. 



(10) 



Therefore, resonances occur when B kk , =^ and the total 
proper time T takes discrete values 



T„ 



2mr 



(11) 



LU k + UJ k > 

We emphasize that the value of T n does not depend 
the segment, i.e. &n = <rf . After simple calculations it on the details of the travel scenario; however, the total 
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amount of entanglement generated EV = NB^, does 
depend on the specifics of the trajectory thorough the 
Bogolubov coefficient. It is also possible to calculate the 
average number of excitations (N k ) in each cavity mode 
while on resonance. We find with the result is also pro- 
portional to 

Sample travel scenario. — We now specialize to our 
sample travel scenario which corresponds to a cavity 
which is initially inertial, travels with proper accelera- 
tion a = h/L for proper time r, coasts for proper time 
t, travels with proper acceleration (in either the same or 
opposite direction) a' = h'/L for proper time r and fi- 
nally coasts for proper time t. The total proper time of 
the segment is T n = 2(r + t) and the symplectic trans- 
formation is 



S kk ' =U kk > (t)V kk \ {ti)U kk , (r)V kk , (ti) x 
U kk ,{t)V^,{h)U kk ,{T)V kk ,{h). 



(12) 



In this case, we find an analytical expression for the first 
order correction to the beta Bogolubov coefficient, 

\B$>, | - c kk , 1 1 - g* k gt, ||1 + eyglgl, ft ft, \h (13) 

where the phases are given by g k — exp(zwfcr) and 
f k = exp(iui k t) and we chose a' = eya (y > 0). The 
parameter e is equal to 1 when the acceleration of the sec- 
ond period has same the direction as a and —1 when the 

direction of acceleration changes. c kk > = ~^^5^p p]3 — ~ 
are coefficients which vanish for evenly separated modes. 



Substituting Eq. (Ill in Eq. (13) we find that the loga- 



rithmic negativity at resonant times is, 

= NC kk ,\(l-(-l) n e^ + ^l + (-irey)\h. (14) 

Note that the logarithmic negativity vanishes when n is 
even and the time of coasting is t = 2i:m/ (ui k + oj' k ) and 
when n is odd and t = (2m + l)ir/{oj k + u> k ) with m G N. 
In the case the accelerations have the same magnitude 
(y = 1) maximal amounts of entanglement are generated 
when n is even and t — n(2m + l)/(u) k + u' k ) in the case 
when the accelerations have the same direction (e = 1) 
and when n is odd and t = 2-Km/(oj k + u>' k ) when the 
accelerations alternate direction (e = —1). This behavior 
is shown in Fig.[T]where we plot the symplectic eigenvalue 
pjy after N — 5 segment repetitions as a function of the 
proper time of acceleration r and the time of coasting t. 
We considered a cavity length of L = 1, massless modes 
k = 1 and k' — 2 and accelerations a = i 
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Interestingly, the special case of alternating accelera- 
tion directions (e = —1) and no coasting (t = 0) corre- 
sponds to the standard dynamical Casimir setting where 
the cavity oscillates periodically as a whole. A resonant 
enhancement of particle creation occurs in the dynami- 
cal Casimir effect [llj which was recently demonstrated 
in the laboratory in a superconducting circuit consisting 




FIG. 1. The correction to the symplectic eigenvalue pffi after 
N — 5 segment repetitions as a function of the proper time 
of acceleration r and the time of coasting t. We considered a 
cavity length of L = 1, massless modes k = 1 and k' — 2 and 
accelerations a = a' = 10 -4 . 



of a coplanar transmission line with a tunable electrical 
length which produces an effective moving boundary [TI5] . 

Discussion. — We have shown that non-uniform mo- 
tion can generate two-mode quantum gates which pro- 
duce observable amounts of entanglement. Finding ways 
to create significant amounts of entanglement in rela- 
tivistic settings is of great interest since entanglement is 
necessary for quantum communications and information 
processing [20 . Recent studies in relativistic quantum 
information show that small amounts of mode entangle- 
ment are created when a cavity undergoes non-uniform 
motion. [S] . We show that particle creation and bi- 
partite mode entanglement can be linearly enhanced by 
repeating any travel segment periodically. Via the equiv- 
alence principle our results suggest that fluctuations of a 
gravitational field can produce entanglement. For exam- 
ple, consider a small cavity containing a bosonic field in 
its vacuum state free falling in the presence of a gravi- 
tational field [21]. Entanglement between the modes is 
generated by suddenly holding the cavity at a fixed po- 
sition against the action of the gravitational field. If the 
cavity's position changes periodically or the gravitational 
field fluctuates, the entanglement can be enhanced. 
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